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ABSTRACT 

Debris disks with extremely large infrared excesses (fractional luminosities > 
1CT 2 ) are rare. Those with ages between 30 and 130 Myr are of interest because 
their evolution has progressed well beyond that of protoplanetary disks (which dissi- 
pate with a timescale of order 3 Myr), yet they represent a period when dynamical 
models suggest that terrestrial planet building may still be progressing through large, 
violent collisions that could yield large amounts of debris and large infrared excesses. 
For example, our Moon was formed through a violent collision of two large proto- 
planets during this age range. We report two disks around the solar-like stars ID8 and 
HD 23514 in this age range where the 24 /im infrared excesses vary on timescales of 
a few years, even though the stars are not variable in the optical. Variations this rapid 
are difficult to understand if the debris is produced by colli sional cascades, as it is for 
most debris disks. It is possible that the debris in these two systems arises in part from 
condensates from silicate-rich vapor produced in a series of violent collisions among 
relatively large bodies. If their evolution is rapid, the rate of detection of extreme 
excesses would indicate that major collisions may be relatively common in this age 
range. 

Subject headings: circumstellar matter — infrared: stars — stars: individual (2MASS 
J08090250-4858172, HD 23514) — planets and satellites: formation 



1. Introduction 

Circumstellar disks are key signposts for planetary system evolution both because they are 
associated closely with planets, and because they are much more readily observed than the planets 
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themselves. Opticall y thick proto planetary disks are born around young stars and dissipate over a 



time < 10 Myr (e.g. JWyatt 



2008). Afterwards, optically t hin planetary debris disks can emerge, 



sustained by the fragmentation of colliding planetesimals (|Wyattll2008|) . Debris disks are readily 
detected through their infrared excess emission (over the stellar photospheric output), which is 
produced by dust warmed by the star. Such excesses are found over the entire age range of main 
sequence stars to ages of ~ 10 Gyr. Therefore, they can be used to probe the presence and state 
of planetesimal systems (such as the analogs of asteroid and Kuiper Belts) and, indirectly, the 
evolution of the accompanying planetary systems. In particular, debris disks can be used to search 
for phases in the evolution of other planetary systems that played major roles in the evolution of 
the Solar System and in the formation of terrestrial planets. 

One such phase extends roughly from 30 to 130 Myr, well separated from the era o f proto- 
planetary disks but when dynamic al models (e.g..lChambersll200ll: iMorishima et al.ll2010 ) predict 
that giant impacts will still occur (Right er & O'Brienll201l|) . For example, our Moon formed dur- 
ing this perio d through a massive collision between the proto-Earth and a Mars-sized proto-planet 
(ICanupll2004|) . Such events should yield huge amounts of debris; indeed, a few candidates have 
been found. They take the form of infrared excesses from heated debris that far exceed th ose ex- 
pected from the trend for quiescently evolving systems (see summary in lBalog et al.ll2009|) . These 
extreme systems (defined as total emission > four ti mes the stellar ph otospheric output at 24 /mi) 
occur around only ~ 1% of solar-like (F5-G9) stars (|Balog et al. 2009 ). If evolutiona ry timescales 



simi lar to those of conventional debris systems are adopted (on the order of 3 Myr, iGrogan et al 



200 1|) . then the low level of incidence of such excesses implies that major collisions are rare. 



In this paper, we present Spitzer data that show that two of these extreme systems, ID8 in NGC 
2547 and HD 23514 in the Pleiades (both with fractional debris disk luminositie s > 2 x 10~ 2 ), vary 
signi ficantly at 24 /im over a few years. The age of NGC 2547 is ~35 Myr (|Jeffries & Oliviera 
2005 ). while that of the Pleiades is ~ 1 20 Myr ( Stauffer et al.lll998 ). Despite their large and variable 
infrared excesses, both ID8 and HD 23514 are well past the age of protoplanetary disks and the 
variations we report ar e unlikely to be related to those found to b e ubiquitous in the latter disk type 
at ages < 10 Myr (e.g. JEspaillat et al.ll201 lUFlaherty et al.ll2012h . iMelis et al.l d2012h have reported 
a third example of strong variability in an extreme infrared excess in the same 30 — 130 Myr age 
range, so the phenomenon appears to be characteristic of this particular class of circumstellar disk. 

ID8 is a type G6-G7 dwarf (N. Gorlova et al., in preparation) at a distance of ~450_p_c_and 



with optical colors indicating negligible extinction. HD 23514 is of F5V spectral type (|Gray et al 
200 lb at a distance of ~130 pc, with norm al colors for this type (M endozalll967l) when corrected 



for the foreground reddening (jCernislll987|) . ID8 is possibly a binary star based on a radial velocity 
measurement 8 km s _1 different from the cluster mean (N. Gorlova et al., in preparation); a late-M 
companion to HD 23514 has recently been discovered at a projected distance of 360 AU from the 
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star (R odriguez et al l2012h . Neither star is variable in the optical, to within the ~ 1% measurement 



errors. 



2. Data 



Spitzer/MlPS observations of ID8 at 24 /im were obtained in programs for G. Rieke (PID 58) 
and C. Lada (PID 20124), The data for HD 23514 are from programs led by M. Meyer (PID 148) 
and G. Rieke (PIDs 30503, 30566). We also show in Figure 1 IRS spectra; that for ID8 is from 
PID 40227, PI G. Rieke, while that for HD 23514 is from PID 50228, PI I. Song. In summary, we 
have three epochs of MIPS 24 /im observations for both targets, supplemented with a 4th epoch 
observation using IRS spectra. In addition, one epoch of IRAC observation is available for each 
target, from PID 58 for ID8 (PI G. Rieke) and PID 50228 for HD 23514 (PI I. Song). 



The 24 jim data were reduced with the MIPS instrument team Data Analysis Tool ( Gordon et al 



20051) . In addition, a second flat field const ructed from the 24 jim d ata itself was applied to remove 
scattered-light gradients and dark latency (|Engelbracht et al.ll2007l) . We extracted the photometry 



using PSF (point spread function) fitting. The input PSFs were constructed using observed cali- 
bration stars and smoothed STinyTim model PSFs jj, and have been tested t o insure the photometry 



results are consistent with the MIPS calibration (|Engelbracht et al. 2007 ). Th e final photometry 



error s also include the errors from the detector repeatability (< 1% at 24 /im, lEngelbracht et al. 



2007|) . These errors sho uld apply directly to HP 23514, which is similar in brightness to the stan- 
dard stars described in lEngelbracht et al.l (|2007|) . We measured an isolated, similar brightness star 
~4' south of ID8 to test the 24 fim photometry stability among the three epochs of data. The 
measured 24 fim flux (5.46 mJy) is constant within 1.5% (rms), consistent with the repeatability 
of MIPS photometry on brighter stars. The measured flux densities and 24 fim magnitudes ([24]) 
assume 7.17 Jy as the zero magnitude flux density. 



The IRAC photometry of ID8 has been published by iGorlova et all (120071) . For HD 23514, 
the IRAC BCD images were retrieved from the Spitzer archive, and mosaicked with the MOPEX 
software with a scale of 0.6" pixel -1 . The final mosaic images were checked against the archival 
PBCD images for consistency. Aperture photometry was then performed on the verified mosaic 
images with an aperture radius of 15" and sky annulus between 24" and 36" from the opto-center. 
The flux conversion issue in the SI 8. 18 IRAC processing^ is corrected. 



IRS spectra were reduced and extracted using the SMART software (Hi gdon et al.l 12004; 



'Krist, J., 2006, Tiny Tim / Spitzer User's Guide, Version 2.0 
2 http://irsa.ipac.caltech.edu/data/SPITZER/docs/irac/sl8.18fluxconv.shtml 
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Lebouteiller et al.ll2010|) . We also computed synthesized 8.0 and 24 /mi photometry by integrating 
the IRS spectra and comparing with similar integrals over an A-star spectrum. We adopt 5% uncer- 
tainty for the sy nthesized photometry which is co nsistent with the cross calibration among IRAC, 
MIPS, and IRS iCushing et al. 2006 : Carey 2010 ). The 24 /mi photometry is summarized in Table 
1 and shown in Figure 1. Both stars are found to vary at 24 /mi at a high degree of significance. 

Visible photometry was collected from the literature (Table 2). Neither star has been noted 
to vary; in fact, HD 23514 has been used as a local photometric standard to study variations 
in neighboring stars (jAlphenaar & van Leeuwenl Il98l|) . Because there were only four relevant 
measurements for ID8, we obtained two more; eight measurements were already available for HD 
23514. The data for both stars show a scatter of only ±1%, consistent with measurement errors 
and corroborating the lack of variability at visible wavelengths. 

WISE observations of ID8 were performed in mid May, 2010, nearly 6.5 years after the IRAC 
observation. Here we consider only the Wl and W2 bands because of their spectral proximity to 
the two short wavelength bands in IRAC. We randomly select three other NGC 2547 member stars 
whose JHK, 3.6 and 4.5 /im brightnesses are all similar to that of ID8. The (K s - W2) colors 
of the three comparison stars are in a narrow interval between 0.00 and 0.03, while that of ID8 is 
approximately 0.67, confirming the 4.5 /mi excess. On the other hand, the ([4.5] - W2) colors of 
the comparison stars range from 0.01 to 0.05, while that of ID8 is 0.28. Given the similarity of the 
IRAC 4.5 /mi and WISE W2 bands, the significant difference between the fluxes in these bands 
at different epochs is most easily explained as the result of variability. The same trend may also 
be observed at 3.6 /im (or Wl) at a low degree of significance. However, variations in the 3-4 
/mi region are not seen for HD 23514 over the ~1.5 year time span between the IRAC and WISE 
observations. 

In summary, we found 10% — 30% peak-to-peak variation at 24 /im on yearly timescales 
while no changes were found at optical wavelengths for these two systems. The change between 
IRAC and WISE, as well as the IRS synthetic photometry, suggests that variations are also possibly 
seen at 4.5 and 8.0 /mi. 



3. Discussion 



HD 23514 and a few other of extreme-excess stars (e.g., HD 15407. [Fujiwara et al.l 12009 : 
Melis et al.llioioh show s pectra in the 10 / mi region (see Figure 1) indic ative of the presence o f 
finely divided silica dust (fLisse et al.ll2009|) or analogs of silica "smoke" (|Kimura & Nuthll2007|) . 
The spectrum of ID8 is also complex, with features indicating significant amounts of finely di- 
vided crystalline silicates (N. Gorlova et al., in preparation). The large equivalent widths of the 
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features in these spectra require grain sizes of order 0. 1 /mi (IBouwman et al. 2008), well below the 
blowout size for these stars. The dust p roperties suggest t hat the infrared excesses of these stars are 



associated with violent, recent events (|Lisse et al. 



2009) and hence may differ signi ficantly from 



typical planetary debris systems whose mid-infrared spectra tend to be featureless (Wyatt 2008) 



and whose evolu tion is dominat ed by an overall slow monotonic decay with time (|Rieke et al. 



200alWyattll2008h . BD +20 307 (ISong et alJl2005L several Gyr old), which also has a very strong 
emission band around 10 /mi, may also vary by ~7% in the WISE W3 and WA bands (12 and 22 
/im) between two scans ~188 days apart, although it is probably stable at Wl, suggesting that the 
mid-infrared variability might be correlated with the presence very fine dust particles. 

Traditional planetary debris disks are sustained by collisional cascades in which populations 
of planetesimals collide until they are ground down to micron-sized dust, when they are ejected 
by radiation pressure or spiral into the star due to Poynting -Robertson drag. The timescale fo r the 
overa ll changes in the infrared excesses in these systems is expected to be mill ions of years dwvatt 
2008). However, around stars of solar mass or less, theoretical models by Kenyon & Bromley 
(|2005l) indicate that short-term spikes in the 24 /im output can occur due to large collisions (of 
~100 km bodies). Thus, the variations support the arguments that we are seeing the consequences 
of individual collisional events and the resulting rapid generation of large swarms of particles. 



We can estimate the timescales for such events from the scaling relations given by IWyatt 
d20Q8h . To apply these relations, we need to estimate the location of the planetesimal belts around 
the two stars. Because small grains associated with recent transient events dominate the mid- 
infrared properties, it is not possible to derive the properties of the parent body planetesimal sys- 
tem from the mid-IR measurements. Instead, we assume that the parent body populations are 
similar to those around stars of similar spectral type and age. The infrared excesses at 5.8 and 
8.0 /mi trac e dust within ~1 AU of a solar-type s tar. W e have complemented studies of these 
exces ses by iGorloya et al] (|2007|) : ICarpenter et al.l (|2009l) by analyzing the data for the Pleiades 



from Stauffer et al. 



(2007). In the latter case, we found that none of the stars with [3.6] < 10 
(corresponding to solar-like absolute magnitudes) have excesses. The net result is that, in the age 
range (30 - 130 Myr), there are only two cases other than ID8 and HD 23514 with reliably detected 
excesses at 5.8 /im and 8.0 /mi out of a total of 352 solar-like stars observed. This result shows 
that planetesimal belts in d ebris systems in this age range are nearly always located outside 1 AU. 
If into equations 14 - 16 of lWyattl (I2008|) . we substitute solar masses and luminosities, a fractional 
luminosity of 2 x 10~ 2 , a radius of 1 AU, dr/r of 0.5, and eccentricity of 0.2 (values to yield as 
rapid evolution as possible), we find that the time to remove a 100 km diameter body by collisional 
cascade is t c > 100 years. On the other hand, the time to remove 1 mm particles is t c ~ 0.1 year, 
an order of magnitude shorter than the time span of the observations. This suggests that small par- 
ticles are replenished by secondary collisions. Given that the variations are ~10% of the fractional 
luminosity, the mass associated with them is ~10~ 2 M®, integrating the particle size distribution 
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from 0.1 /im to 100 km. Significantly lower mass estimates are only possible if the planetesimal 
belt is well inside 1 AU. 



Grogan et al.l (|200ll Figure 19) present detailed models that find a characteristic timescale of 
3 x 10 6 years for the decay of dust from collisions in the asteroid belt. Th e behavior is supporte d 



by the role of recent collisions in the production of zodiacal dust bands (|Nesvorny et al. 



2003) 



This result provides an independent test of the timescale for decay of excess due to collisional 
cascades around ID8 and HD 23514. For the population of disks around stars of s imilar mass and 
age to these two, a typical fractional luminosity in a quiescent state is 5 x 10~ 4 (|Carpenter et al 



2009J). Most of this emission is from cool dust that must be well outsid e a few AU from the star. 
The fractional luminosity for warm dust is an o rder of magni tude lower ( Carpenter et al.ll2009 ). i.e. 
5 x 10~ 5 . These values, the scaling relations in IWyattl (|2008f) . a solar fractional luminosity of 1CT 7 
(IBackman & ParesceHl993|) . and orbital radius of 2.5 AU for the zodiacal cloud and asteroid belt 
indicate ~100 years as the characteristic time for decay of a transient event via collisional cascade 
at 1 AU around ID8 or HD 23514. This value agrees well with that obtained directly from the 
scaling relations. 

Although this estimate is rough, it indicates that generating the variable excesses through col- 
lisional cascades may not be impossible, but would require extreme assumptions. We therefore 
will explore alternatives. The possibility that the 24 fim changes reflect varying rates of heating by 
the star is made unlikely by the complete lack of variations in the visible for both stars. Phenomena 
that avoid any kind of collisional cascade will inherently have faster timescales. One such process 
might be a collision that caused a large body to lose its regolith, injecting pre-processed small dust 
into the interplanetary space around the star. However, very small grains tend to be a bsent in plan 



etesi mal regoliths because radiation pressure slowly drives the small particles away (|Masiero et al 



20091) : the regolith hypothesis is possibly contradicted by the evidenc e for small grains in the mid - 
IR spectra. In addition, regoliths have a full size spectrum of objects (IHousen & Wilkeninglll982|) . 
so such a regolith ejection event would be analogous to launching a fully developed collisional 
cascade and it would have a similar evolutionary timescale to the cascades already considered. An 
alternative might be the disintegration of one or several extremely large comets. The most finely 
divided dust released would have a short dwell time near the star, since it would be ejected by radi- 
ation pressure or lost through Poynting-Robertson drag; the larger bodies would remain for some 
time because the comet orbit would not coincide with a dense planetesimal belt and would have 
inadequate density of objects for vigorous collisional activity. The timescale for the first process is 
probably too short, making the resulting excesses too rare to agree with the simultaneous detection 
of a number of them, and that for the second too long to explain the variability of ID8 and HD 
23514. 



None of these possibilities can convincingly be shown to produce the huge amounts of dust 
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around these stars while having the necessary short timescales for variability. Therefore, we con- 
sider another possibility - that the dust condenses from silica-rich vapor associated with violent 
collisions among large bodies. For example, in massive collision s such as the one that formed our 
Moon, of order 20% of the mass emerges as vapor rich in silica (|Canupll2004|) . The production of 
silica-rich vapor can occur at collisional velocities above about 2 km s _1 based on simple energetic 
arguments. However, a threshold of about 10 km s 1 is generally ado pted because a substa ntial por- 
tion of the energy is deposited in collisionally produced fragments (IHornung et alj|2000h. Labora- 
tory e xperiments suggest that the vapor will condense quickly into silica "smoke" (IKimura & Nuth 



20071) . and forsterit e-like crystals are the most common crystalline form to grow from the con- 



densation products (|Kobatake et al.ll2008|) . The spe ctrum of silica sm oke (|Kimura & Nuthll2007r) 
resembles that of the infrared excess of HD 23514 ( Rhee et al.l l2008). while the small grains that 
dominate the spectrum of ID8 are rich in forsterite-like materials (N. Gorlova et al., in prepara- 
tion). The resulting small particles will be cleared quickly by radiation pressure and the Poynting- 
Robertson effect. The cooling and di sk spreading timescale for the silica-rich pro tolunar disk that 
later formed the Moon is ~one year (Thompson & Stevensonlll988l: ICanupl l2004). However, most 
of the ejected mass in the moon-forming impact remained grav nationally bound to the Earth; the 
small solid angle viewed from the Sun indicates that such a single event could not raise the infrared 
luminosity high enough to explain the excesses in these systems. However, the significant amount 
of escaping material from a major impact might initiate a violent collisional cascade and elevate 
the infrared excess of the star for a longer period. Thus, the extreme excesses around ID8 and HD 
23514 are unlikely to be due to an individual violent collision, but to reflect extended episodes of 
such collisions possibly triggered by a single event. 



4. Conclusion 

Dynamical models indicate that giant impacts continue to build terrestrial planets into the 30 
- 130 Myr age range. A small number of stars have extreme infrared excesses and mid-infrared 
spectra that suggest that such impacts have occurred recently around them. However, the duration 
of these indicators, and hence the frequency of events required to account for the numbers we 
detect, is not clear. We report that the excess is variable on year-timescales around two of these 
stars , ID8 in NGC 25 47 and HD 23514 in the Pleiades; a third variable example has been reported 



by iMelis et al.l (I2012|) . The variability is too rapid to be explained readily if the excesses are 
generated in collisional cascades, as is the case for most planetary debris disks. Instead, it is 
possible that the dust producing the excess condenses from vapor generated in violent collisions. 
The timescale for variations of years, rather than millions of years, suggests that the evolution of 
the infrared excesses might be rapid. In that case, major collisions would need to occur frequently 
to account for the number of stars in this age range that are observed to have extreme excesses. 
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Table 1 . Infrared Photometry by Spitzer and WISE 



Date 


3.6 fim (mJy) 


4.5 fim (mJy) 


5.8 fim (mJy) 


8.0 fim (mJy) 


24 fim (mJy) 


source 


ID8 


2003 Dec 4 

2004 Jan 29 
2006 May 9 

2006 Dec 4 

2007 Jun 16 
2010 May 13-21 


9.74 ±0.1 
12.05 ± 0.26 b 


7.35 ±0.1 
9.11 ±0.17 b 


6.56 ±0.1 


7.64 ±0.1 
6.38 a 


6.23 ±0.10 
8.34 ±0.12 
7.32 ±0.11 
6.74 a 


Gorlova et al. (2007) 
AOR 43 18976 
AOR 16798464 
AOR 16800512 
AOR 21755136 
WISE All-Sky Catalog 


HD 23514 


2004 Sep 22 
2007 Feb 25 

2007 Sep 17 

2008 Sep 18 
2008 Oct 1 

2010 Feb 12-13 


206.3 ±1.1 
225.3 ±4.5 b 


180.7 ± 1.2 
180.6 ±3.1 b 


154.6 ± 0.7 


203.9 ± 1.0 
169.3 a 


67.7 ± 0.7 
68.0 ± 0.7 
61.6 ±0.7 

64.8 a 


AOR 5320192 
AOR 18303232 
AOR 17656576 
AOR 25790208 
AOR 25789952 
WISE All-Sky Catalog 



a Synthetic photometry based on IRS spectrum and the resp onse curves of IRAC and MIPS. The cross-calibration uncertainties of 
the Spitzer instruments are expected to be ~ 5% ( Carey 2010). 

b Flux densities at WISE Wl and W2. The effective wavelengths are 3.35 and 4.60 fim, slightly different from those of IRAC. 
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Table 2. Visible 
Photometry 



V I a references 



ID8 



13.14 


12.34 


1 




12.34 


2 


13.12 


12.34 


3 




12.32 


4 


13.12 




5 


13.14 




5 


HD 23514 


9.42 




6 


9.42 


8.67 


7 


9.44 




8 


9.43 




9 


9.42 




10 


9.44 




11 


9.43 




12 


9.44 




13 



a Cousins filter for ID8 and Johnson filter 
for HD 23514. 



References. — (1) Naylor et al. 2002; (2) 
Jeffries et al. 2004; (3) Lyra et al. 2006; (4) 
DENIS catalog; (5) this work, observed on 
201 1 Nov 18; (6) Johnson & Mitchell 1958; 
(7) Mendoza 1967; (8) Rufener 1976; (9) 
Alphenaar & van Leeuwen 1981; (10) Cer- 
nis 1987; (11) van Leeuwen et al. 19 86; (12) 
Tycho catalog, transformed from [Bessell 
l2000l) ; (13) Droege 2006. 
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Fig. 1. — Top: IRS low-resolution spectrum of ID8 (smoothed by three pixels), with our MIPS 
photometric points at 24 /im (black) plus a value from synthetic photometry of the spectrum (red). 
A 5% uncertainty is adopted for the synthetic photometry. The response curve of MIPS 24 ^m is 
overplotted for comparison. Bottom: Same, for HD 23514. 



